Ground state properties and dynamics of the bilayer t — J model 
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We present an exact diagonalization study of bilayer clusters of t — J model. Our results indicate a 
crossover between two markedly different regimes which occurs when the ratio J±/J between inter- 
layer and intra-layer exchange constants increases: for small J±/J the data suggest the development 
of 3D antiferromagnetic correlations without appreciable degradation of the intra-layer spin order 
and the d x 2_ y 2 hole pairs within the planes persist. For larger values of J±/J local singlets along 
the inter-layer bonds dominate, leading to an almost complete suppression of the intra-layer spin 
correlation and the breaking of the intra-layer hole pairs. The ground state with two holes in this 
regime has s-like symmetry. The data suggest that the crossover may occur for values of J±/J as 
small as 0.2. We present data for static spin correlation function, magnetic structure factor and spin 
gap. We calculate the momentum distribution and electronic spectral function of the "local singlet 
state" realized for large J±/J, and show that it deviates markedly from that of a single layer, making 
it an implausible candidate for modelling high-temperature superconductors. 
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I. INTRODUCTION 

In connection with the so-called spin-gap phenomenon 
jjj observed in high-temperature superconductors the 
properties of coupled layers of 2D t — J model have 
recently attracted considerable attention In this 

work we present results obtained by exact diagonaliza- 
tion of small clusters, which may provide a rough guide- 
line for this problem. Our data were obtained by diago- 
nalization of clusters of size 2x8 and 2 x 10, essentially 
the limit that can be reached by the exact diagonaliza- 
tion technique. Clearly these are still very small sys- 
tems, and subtle low-energy effects certainly will escape 
our study. However, one may expect that even small 
clusters can mimick the short range correlation functions 
(which solely determine the total energy) roughly correct, 
and when interpreted with care our calculations thus may 
provide a crude phase diagram of the bilayer. Our key 
result is that already for fairly moderate values of the 
inter-layer exchange a qualitatively new type of ground 
state is realized, which deviates significantly from that of 
a single layer. 

The model under consideration reads 

<i,j>,<7 <i,j> 



H = 



Here the Si are the electronic spin operators, 
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and the sum over < i,j > stands for a 



summation over all pairs of nearest neighbors. Our sys- 
tem consists of two planes, labelled A and B for later ref- 
erence, the z-axis of the coordinate system is taken to be 
perpendicular to the planes. We distinguish between the 
hopping integral and exchange constant between nearest 
neighbors in the same plane, which we denote by t and </, 
and between nearest neighbors in different planes, which 
we denote by t± and J±. There is presently an unclear 
situation as to the correct choice of parameter values, ex- 
perimental estimates for the ratio Jj_/ J vary from 0.085 
over 0.3 Q, up to 0.55 ||. As for the inter-layer hop- 
ping integral t±, the situation is even more unclear. Let 
us note that inter-layer hopping integrals derived from 
band structure calculations for cuprate superconductors 
are meaningless on the level of the t — J model, because 
the "hole" in the t — J model stands for an extended 
object, the Zhang-Rice singlet, which has very different 
hybridization matrix elements with orbitals outside the 
plane than a "bare" electron ||. We will therefore con- 
sider different ratios of the various parameters and study 
their impact on e.g. the nature of the ground state. 



II. UNDOPED BILAYER 

To begin with, we briefly discuss the spin correla- 
tions of an undoped bilayer. Figure |l| shows the nearest- 
neighbor static spin correlation function in plane and 
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out-of plane as well as the spin structure factor for mo- 
mentum tranfer (tt, 7t,0) and (7r, 7r, tt). As J±/J in- 
creases, there is initially a fairly symmetric splitting be- 
tween S(tt, 7t, 0) and S(w, tt, n), their average S(tt,tt) — 
2(5(77, 7T, 0) + 5(7r, 7T, 7r)) remaining almost constant. The 
latter quantity measures the antiferromagnetic correla- 
tions within a single plane. In this range of J±/J, there 
is virtually no degradation of the antiferromagnetic order 
within the planes. The maximum of S(tt, tt, tt) is reached 
rather precisely for J±/ J = 1, and for this value the near- 
est neighbor spin correlation function within and between 
the planes are equal. Quite obviously, for this parame- 
ter value the bilayer is closest to a 3D Heisenberg an- 
tiferromagnet. When J±/J is increased further, S(tt, tt) 
decreases and so does the intra-layer nearest neighbor 
spin correlation function. The spin correlations within 
the plane now decrease and local singlets along the inter- 
layer bonds start to dominate. However, the intra-plane 
correlations still are relatively " stiff' , so that large values 
of J±/J are required to bring them close to zero. 



III. DOPED CASE 

We next consider the doped bilayer, more precisely the 
ground states with two holes of the 2x8 and 2 x 10 
bilayer. These show the same overall trends as the un- 
doped bilayer, however within a much smaller range of 
J±/J. Figure |^ shows the nearest neighbor spin corre- 
lation functions as functions of Jj_/J and demonstrates 
that the crossover to a state with almost vanishing in- 
plane spin correlations and dominant z-axis spin corre- 
lations is much more abrupt than for the undoped sys- 
tem. Hole doping obviously makes it easier for the inter- 
plane terms to dominate over the in-plane exchange. The 
strong difference as compared to the undoped case more- 
over shows that results for undoped bilayer systems may 
have limited significance for the doped bilayer. 
The second major difference as compared to the undoped 
bilayer is, that the crossover from the small to large J±/J 
regime is now accompanied by a ground state level cross- 
ing. For small values of J±/J the 2 hole ground state of 
the 2 x 10 bilayer belongs to the B\ (or d x i_ y i) represen- 
tation of the Ci v point group , reflecting the well-known 
fact that two holes in the t — J model always form 
a bound state with this symmetry. In this state, the two 
holes can be found in the same plane with high proba- 
bility i.e. the in-plane bound state obviously persists to 
some degree. This changes completely when J±/J is in- 
creased beyond a certain critical value: the ground state 
now has the A\ (or s) symmetry and the two holes are 
in different planes almost with probability 1. A rough 
"phase diagram" of the 2 x 10 bilayer is given in Fig- 
ure Table Rjshows the hole density correlation func- 
tion g(R) = 2^j((l — n.i)(l — n i+ g)) for both B\ and 
A\ state near the crossover. An analogous level crossing 
also can be found in the 2x8 cluster, for comparable val- 



ues of J±/J, so that it does not seem to originate from 
some spurious effect due to the special geometry of one 
of these two clusters. Strong inter-layer spin correlations 
thus are incompatible with c£-wave pairing of holes within 
the planes, and change the symmetry properties of the 
ground state. It is interesting to note that the incom- 
patibility of d-wave pairing within the planes and strong 
interlayer spin correlations has also been found by Liecht- 
enstein et al. [jllj within a weak coupling framework, as 
well as by Normand et al. j^] within slave-boson mean 
field theory. This seems to be a very general property of 
bilayer systems. 

At this point, the possibility of fairly strong finite-size ef- 
fects on the critical value of J±/ J becomes obvious: since 
the transition to the state with pronounced z-axis spin 
correlations implies the breaking of the in-plane pair, it 
follows that an overestimation of the in-plane hole bind- 
ing energy will lead to an artificial stabilization of the 
"conventional" B\ ground state. To see this, we first de- 
fine the binding energy as Ag = {E^ h ^ — E^ 01 ^) — 2 • 
(E (m) _ E (v h)) = {E (2 h) + E (v h)) 2 _ E (i h} where E{ uh) 

denotes the ground state energy with v holes. In small 
clusters of t — J model this quantity is negative |l(| , so 
that in the limit J±,t± — > the 2 hole ground state of 
the bilayer will approach a product of an undoped layer 
and a layer with two bound holes. The binding energy 
Ab then represents a kind of charge transfer energy, i.e. 
it gives the cost in energy upon transferring one of the 
two holes to the other layer. Hole binding within the 
planes thus opposes the formation of the A\ state, and 
since it is known that finite-size calculations for the ID 
t — J model tend to substantially overestimate Ab [[l2| , 
this may have an impact on the crossover between the 
two different ground states. In order to further investi- 
gate this issue, we add to the Hamiltonian a repulsive 
interaction between holes in the same plane: 

i.jeA i,j£B 

Such a term does not affect the intra-plane dynamics; 
it only reduces the loss of energy which occurs when a 
hole is transferred between the two layers. In particular, 
if we adjust V — Ab this term precisely cancels the ef- 
fects of the in-plane binding energy which may contain 
large finite-size effects. Figure ^ then shows (for the value 
J/t = 0.4) the change of the crossover value of ( J± / J) cr it 
when the in-plane repulsion is increased. There is an ob- 
vious linear dependence and for completely "balanced" 
intra-layer binding energy (i.e. for vanishing hole trans- 
fer energy) we find (Jj_/J) C rit — 0.17. For J/t = 0.2 the 
same procedure gives (J±/J) C rtt = 0.15 (for completely 
balanced binding energy), so that in the parameter re- 
gion of possible relevance for cuprate superconductors 
the crossover to the unpaired state may occur for fairly 
small values of Jj_/J. 

For completeness we have to address possible ground 
state degeneracies. For the single 10-site cluster with 
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periodic boundary conditions there exists a special per- 
mutation of the sites which leaves all nearest neighbor 
relations invariant; this is analogous to the well-known 
2 hypercube symmetry of the 4x4 cluster with periodic 
boundary conditions JlO| | . This artificial symmetry oper- 
ation therefore commutes with any Hamiltonian involv- 
ing only nearest neighbor terms and introduces artificial 
degeneracies. Clearly, for the 2x10 bilayer there exists an 
analogous permutation of z-axis bonds, leading to artifi- 
cial degeneracies as well. We have always neglected the 
degenerate ground states with finite momentum, because 
we believe that the ground state with vanishing momen- 
tum is the most natural representative of the infinite sys- 
tem. In the single layer case, this is confirmed by the fact 
that the k — ground states of different clusters all have 
very analogous properties, whereas the degenerate finite 
momentum ground states disappear for larger clusters. It 
is moreover straightforward to see, that these degenerate 
states cannot add much new information: for example 
the average nearest neighbor spin correlation functions 
can be expressed as derivatives of the total energy with 
respect to the exchange constants J and Jx, the proba- 
bility to find the two holes in different planes is related 
to the derivative with respect to the parameter V in (Q) , 
so that these quantities must be rigorously the same for 
all symmetry-degenerate ground states. 
We now turn to the spin correlations, shown in Figures 
|5| and |(]. The B\ state shows the same symmetric split- 
ting of S(tt, 7r, 0) and S(tt,tt,tt) which was seen already 
in the undoped spin-bilayer; the inter-plane spin correla- 
tions increase only moderately. Quite obviously, in this 
state the development of inter-layer spin correlations is 
accomplished by the " lock in" of the directions of the lo- 
cal staggered magnetization Ms in the two planes, with 
essentially unchanged intra-plane correlations. By con- 
trast, the Ai state shows dominant inter-plane correla- 
tions and, as shown in Figure |^, already for values of 
J_l/J ~ 1 it rapidly approaches the limit of a product 
state of z-axis singlets. 

In order to further clarify the mechanism by which inter- 
layer spin correlations are built up, we now study the 
probability distribution of the staggered magnetization 
in the ground state. More precisely, we denote by 
the set of basis states which have staggered magnetiza- 
tion n in the A-plane and staggered magnetization v in 
the B-plane. Then, we can define the operator 

r^= l*«X*a|, 

and its ground state expectation value S(fi,u) = (V^u). 
It may be interpreted as the total weight of states with 
single-layer staggered magnetizations fi and v. Figure |?j 
shows this quantity for two values of Jj_/J. For vanish- 
ing inter-layer coupling the ground state is simply the 
product of two single-layer ground states, and if one 
chooses fixed z-component of the momentum, S(/i, v) 
must be completely symmetric, i.e. S(fi, v) = S(fi, —v) = 



S(—[i, v) = S(— /i, —v). The data then show, that for in- 
creasing Jj_/J there is first a slight shift of weight from 
the (+, +) and (— , — ) quadrants to the (+, — ) and (— , +) 
quadrants: the staggered magnetization of the two layers 
"locks in". However, states with large staggered magne- 
tization still do have appreciable weight, and the correla- 
tion between the values of Ms in the two planes is not yet 
strong. This changes completely when we increase J±/J 
beyond the crossover to the A\ ground state: S(fi, v) is 
now concentrated near the line fi = —v, and along this 
line it drops sharply towards large values of Ms- The 
values of Ms in the two layers are now strongly corre- 
lated: when there is e.g. a fluctuation towards a large 
positive value of Ms in the A-plane, it must be accom- 
panied by a fluctuation to a large negative value in the 
.B-plane. These results are precisely what one would ex- 
pect if the spins were distributed randomly within the 
planes, but under the constraint that nearest neighbor 
spins in z-direction always be antiparallel. Analysis of 
the probability distribution S^^ thus confirms the sce- 
nario inferred from the spin correlation function: there is 
a regime of small J±/ J, where the system responds to the 
inter-layer coupling by a weak "lock-in" of the staggered 
magnetization in the two planes, while no appreciable 
change of the intra-layer correlations takes place. For 
large J±/J, on the other hand, we have the clear signa- 
tures of the z-axis singlet state. 
We next consider the " spin gap" , defined as 

A spin = -Eo,r (tt, t, tt) - E 

Here Eq denotes the energy of the ground state (which 
has momentum 0) and Er (tt, 7t, tt) denotes the lowest en- 
ergy in the subspace of states with momentum (tt, tt, tt) 
and the same point group symmetry Tq as the ground 
state. Inspection shows that the k = ground state is al- 
ways a singlet, the (tt, tt, tt) state always a triplet, so that 
A sp i„ indeed represents the gap which would be observed 
in the dynamical spin correlation function with momen- 
tum transfer (tt,tt,tt). We note that A sp i n is prone to 
massive finite-size effects, as can be seen e.g. from the 
fact that even an undoped 2D cluster has a "spin gap" 
of w 0.4 J @. The values of ^spin determined from the 
small clusters thus may at best indicate a rough qualita- 
tive trend, but clearly have no quantitative significance. 
Then, Figure ||a shows the spin gap over a wide range of 
Jj_/J. It is remarkable that the gap initially decreases 
with increasing J± and in fact shrinks appreciably at 
the level crossing to the strong z-axis correlation (Ai) 
state. It is only for values of J±/J > 1 that the spin 
gap becomes large and grows linearly with Jj_, reflect- 
ing the saturation of the z-axis spin correlation function 
(see Figure ||) . For an undoped system and in the limit 
Jj_/J — > oo it is easy to see that a single z-axis triplet 
bond can propagate in a " background" of z-axis singlets 
with a nearest neighbor hopping element of J/2, so that 
the spin gap should be J± — z-(J/2), with z the number of 
nearest neighbor z-axis bonds. The resulting asymptotic 
expressions A sp i n = J± — J for the ladder (z = 2) and 
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A spin = J± — 2J for the bilayer (z = 4) indeed provide 
excellent estimates for the undoped systems (compare 
Figure 2 in Ref. ||); for the doped system, however, the 
numerical result rather suggests the asymptotic expres- 
sion A sp i n = Ji — (J/2), which would imply a reduction 
of the effective hopping element for the triplet-bond by 
a factor of 4. This indicates, that the propagation of the 
z-axis triplet is inhibited by the presence of mobile holes, 
so that there must be a strong interplay between the lo- 
cal singlets/triplets and the holes. 

Figure @b shows the spin gap for other parameter values 
in the range J±/J < 1 and demonstrate that increasing 
values of J±/ J initially tend to reduce A spini as would 
be expected from an enhancement of antiferromagnetic 
correlations, ft is interesting to note that the crossover 
B\ — ► A\ never leads to an increase of A sp i n . While the 



significance of these data should not be overestimated, 
we may state that in the small clusters there is no indi- 
cation of any widening of the spin gap due to interlayer 
coupling, at least not for "reasonable" values of J±/J. 



IV. LOCAL SINGLET STATE 

We now proceed to a study of the " z-axis singlet state" , 
which is realized for larger J±/J. As has been shown 
above, this state may appear already for quite moderate 
values of J±/J ~ 0.2, which have in fact been postulated 
to be realized in actual high-temperature superconduc- 
tors. In fact, this state may be considered a qualitatively 
new feature of the bilayer system, as compared to the 
single-layer ground states studied extensively before. We 
first consider the limiting case of vanishing in-plane pa- 
rameters J and t gj. Depending on the relative magni- 
tude of t± and Jj_ , the ground state of the 2xJV bilayer 
with an even number of holes rih can take two quite dif- 
ferent forms (see Figure ^): it is either the product of 
N — n/j/2 z-axis singlets and n/i/2 empty z-axis bonds 
with energy (N — rih/ 2) • J±; or it is a product of N — rih z- 
axis singlets and rih singly occupied z-axis bonds with en- 
ergy (N—nh)-J±—nh-t±. The degeneracy of this ground 
state is lifted by the intra-plane terms of the Hamiltonian. 
We assume that the most important intra-plane term in 
the Hamiltonian is the hopping integral t. Then, a singly 
occupied bond can propagate to a nearest neighbor in a 
single-step process with the effective hopping integral t/2, 
whereas an empty bond propagates to a nearest neighbor 
via a two step process where the intermediate state has 
a broken z-axis singlet; in the regime of large J±/t per- 
turbation theory would give an effective hopping integral 
~ i 2 / Jj_. We thus have the interesting situation that, de- 
pending on the relative strength of parameters, two very 
different "fixed points" may be appropriate to describe 
the bilayer system: if the arrangement of holes in empty 
bonds is preferred, the system should correspond to a gas 
of rih/2 hard core bosons, propagating via nearest neigh- 
bor hopping of strength ~ i 2 /J±. On the other hand, if 



the arrangement of holes along singly occupied bonds is 
preferred, the bilayer should correspond to a system of 
rih spin- 1/2 fermions, propagating via nearest neighbor 
hopping of strength t/2. In both cases, however, there is 
a positive nearest neighbor hopping matrix element for 
either hole- like " Fermion" or " Boson" , so that the holes 
may be expected to accumulate at the in-plane momen- 
tum (tt, tt). 

The calculated values of the electron momentum distri- 
bution function n(k) = (c- Cr ), given in Tables || and 
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then indeed show an enhanced hole-occupancy of 
(tt, 7r) in the A\ state even for moderate values of J± and 
tj_ : this is very pronounced in the 2x8 system, but also 
in the 2 x 10 bilayer one can realize an enhanced depres- 
sion at (tt, 7t) and a tendency towards the "flattening" of 
n(k) in the remainder of the Brillouin zone. This indi- 
cates that the momentum distribution with "hole pock- 
ets" at (tt, tt), which would be expected in the limit of 
large inter-layer coupling indeed seems to persist in the 
A\ state also for moderate values of the inter-layer cou- 
pling. This picture is confirmed by the single particle 
spectral function A(k,u>), shown in Figure |l0|. Due to 
restrictions in computer memory, we have evaluated this 
quantity only for the 2x8 bilayer. The topmost peaks 
in the photoemission spectra are located at (tt/2,tt/2) 
and (tt, 0) (due to the special symmetry of the 2D 8- 
site cluster these momenta are in fact identical for the 
A\ state), but there is no appreciable low energy weight 
in the inverse photemission spectrum at either of these 
momenta: quite obviously these /c-points are " occupied" . 
This is in clear contrast to the spectral function for the 
B\ state |l4j which showed a rather obvious Fermi level 
crossing at (tt,0). Another notable feature is the strong 
contrast between the relatively sharp peaks in the bond- 
ing channel (i.e. k z = 0) and the more diffuse spectra 
in the antibonding one, which have markedly enhanced 
incoherent continua. This indicates that electrons in the 
bonding channel are the "more well-defined" excitations 
of the system, as one would expect it for a product state 
of doubly and singly occupied bonds. The strong differ- 
ence between the k z = and k z = tt spectra is again in 
contrast to their almost identical shape in the B\ ground 
state |Ti| ]. The available data thus rather consistently 
suggest that even for moderate values of J± and t±_ the 
"z-axis singlet state" bears a strong resemblance to the 
J±,t± — > oo limit, which in particular implies the accu- 
mulation of holes at (tt, tt). 



V. SUMMARY AND DISCUSSION 

In summary, we have studied the effects of interlayer 
exchange in a bilayer t — J model. The data rather con- 
sistently suggest a crossover between two very different 
regimes: for small and moderate values of the inter-plane 
exchange coupling there is no qualitative change of the 



4 



spin correlations, the in-plane hole pairing with d x 2_ y 2 
symmetry persists and the inter-plane exchange induces 
3D antiferromagnetic correlations. For increasing values 
of J±/J there is a crossover to a state with dominant 
inter-layer spin correlations, where the intra-layer hole 
pairing is completely suppressed. We note that this fea- 
ture may well persist in the infinite system if the pair- 
ing state remains short ranged; in this case, it would be 
the short range corrrelations which dominate the physics, 
and these may be expected to be described properly even 
in the small clusters. The data finally do not show any 
indication of a widening of the spin gap at the crossover. 
The result that strong inter-layer coupling suppresses the 
in-plane d x 2_ y 2 pairing clearly would be hard to recon- 
cile with the ever growing experimental evidence for d- 
type symmetry of the superconducting order parameter 
in the cuprate superconductors fl5| and the experimen- 
tal fact that T c seems to increase with the number of 
closely coupled layers in these compounds. In addition, 
our data indicate that the " z-axis singlet state" realized 
for large J±/ J has fairly unusual properties such as a ten- 
dency towards " hole pockets" at the in-plane momentum 
(tt, 7r), so that it does not seem to be a very promising 
candidate for modelling real high-temperature supercon- 
ductors. Moreover, the relatively strong difference be- 
tween the z-axis singlet state and the " conventional" sin- 
gle layer ground state would make it hard to understand 
why many properties of single and double layer cuprates 
are relatively similar. 

Clearly, the simplest way to resolve these problems is to 
assume that the values of Jj_/J in the actual materials 
are too small to ever induce the z-axis singlet state. Un- 
fortunately strong finite size effects do not allow us to 
give a very accurate estimate of the critical value where 
the crossover occurs. Then, a possible enhancement of 
antiferromagnetic spin correlations due to the onset of 
3D correlations in the more conventional ground states 
realized for lower Jj_/J may be more easy to reconcile 
with experiment: it seems plausible that pairing theories 
which rely on antifcrromagnetism, such as the antifer- 
romagnetic spin fluctuation theory |l6f| or the antiferro- 
magnetic van-Hove scenario [1^ would profit from the 
onset of 3D correlations, so that the exchange coupling 
of several CuOi planes could be favourable for achieving 
higher T c 's. 
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TABLE I. Comparison of the hole density correlation func- 
tion g(R) in the lowest state with Bi symmetry and the lowest 
state with A\ symmetry of the 2x10 bilayer t — J model. Pa- 
rameter values are t = 1, J = 0.4, t± = 0.5, J±/J = 0.8. The 
value of J± I J is chosen slightly above the crossover Bi — ► Ai , 
so that the Ai state is the ground state, the Bi state being 
0.0289 above the ground state. 





■R-in — plane 


Si 


(0,0) 


(1,0) 


(1,1) 


(2,1) 


R* 
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0.000 


0.0791 
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0.0213 


0.0151 


Ai 


(0,0) 


(1,0) 


(1,1) 


(2,0) 


Rz 
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0.000 


0.0115 


0.0095 


0.0119 


Rz 
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0.2217 


0.1038 


0.0573 


0.0380 



TABLE II. Comparison of the momentum distribution 
function n(k) for the lowest state with Ai symmetry and the 
lowest state with Bi symmetry of the 2 x 10 bilayer t — J 
model. Parameter values are like in Table 1. 
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0.5357 
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0.3132 
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(¥>T?) 
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kz = 





0.5399 


0.5278 


0.4374 


0.3424 


k z = 


7T 


0.5416 


0.4815 


0.3668 


0.3223 



TABLE III. Comparison of the momentum distribution 
function n(k) for the lowest state with Ai symmetry and the 
lowest state with Bi symmetry of the 2x8 bilayer t — J 
model. The value of J±/J = 0.6 is chosen slightly above the 
crossover Bi — > Ai, so that the Ai state is the ground state, 
the Bi state being 0.0154 above the ground state. The other 
parameter values are t = 1, J = 0.4, t± = 0.5. 
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0.5469 
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0.2915 
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7T 
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0.4624 


0.2525 


0.2761 
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FIG. 1. (a) Static nearest neighbor spin correlation func- 
tion {Si ■ Sj) for the ground state of the undoped 2 x 10 bi- 
layer. (b) Static spin structure factors S(n, %, 0) (diamonds) 
and S(7r,7r,7r) (crosses) and average S(n,n) (squares) for the 
same state. 

FIG. 2. Static nearest neighbor spin correlation function 
{Si ■ Sj} for the ground state of the 2 x 10 bilayer with two 
holes. Parameter values are t = 1, J = 0.2, t± = 0.1. 

FIG. 3. "Phase diagram" for the 2 x 10 bilayer with two 
holes for two values of tj_/t. The figure indicates the regions 
of stability of the Ai and Bi ground state for different pa- 
rameter values. The respective critical values of J±/J have 
been obtained by interpolating ground state energies and are 
accurate only to ~ ±0.05. 

FIG. 4. Critical value of (Jj_/ J) c rit for the crossing 
Bi — > Ai as a function of the intra-plane repulsion param- 
eter V. 

FIG. 5. Nearest neighbor spin correlation function {Si ■ Sj) 
for the 2 x 10 bilayer with two holes, (b) Static spin structure 
factors S(n,n,0) (diamonds) and S(n,n,n) (crosses) and av- 
erage S{n, n) (squares) for the same system. Parameter values 
are t = 1, J = 0.4, t± = 0.5. 

FIG. 6. Same as Figure § for J = 0.2. 

FIG. 7. Probability distribution Sp, v for the ground state 
of the 2x8 bilayer with two holes and different Jj_/J. The 
distribution has been averaged over the areas indicated in the 
figure. Parameter values are J±/J = 0.4 (a) and J±/J = 1 
(b), the other parameters are t = 1, J = 0.4, t± — 0.5 for 
both figures. 

FIG. 8. (a) Spin gap Aspin as a function of J±/J for 
J = 0.2, t± = 0.1, 2 x 10 bilayer. 

(b) A 3p i n for different parameter values in the 2 x 10 bilayer, 
t = 1. 
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FIG. 9. Possible ground states of the bilayer model in the 
limit t, J — > 0: the holes are on "singly occupied bonds" (a) 
or on "empty bonds" (b). 



FIG. 10. Single particle spectral function A(k, u) for the 
2x8 bilayer. The vertical line markes the Fermi energy, the 
frquency range uj > Ef {u> < Ef) corresponds to electron 
creation (annihilation). The full (dashed) line refers to the 
bonding (antibonding) channel. Parameter values are like in 
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